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Q"^ , scribing the dynamics of mean fields and their statistical fluctuations in a non- 

^1 \ Abelian plasma in or out of equilibrium. A general procedure is detailed for 

1-^ I integrating-out the fluctuations as to obtain the effective transport equations for 

the mean fields. In this manner, collision integrals for Boltzmann equations are ob- 
^^' tained as correlators of fluctuations. The formalism is applied to a hot non-Abelian 

r~| , plasma close to equilibrium. We integrate-out explicitly the fluctuations with typical 

momenta of the Debye mass, and obtain the collision integral in a leading logarith- 
mic approximation. We also identify a source for stochastic noise. The resulting 



Starting from classical transport theory, we derive a set of covariant equations de- 



\^ • dynamical equations are of the Boltzmann-Langevin type. While our approach is 

based on classical physics, we also give the necessary generalizations to study the 
quantum plasmas. Ultimately, the dynamical equations for soft and ultra-soft flelds 
change only in the value for the Debye mass. 
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I. INTRODUCTION 

This article presents in full detail an approach to the dynamics of non-Abelian plasmas 
based on classical transport theory, the main results of which have been summarized in [0]. 

In the recent years, there has been an increasing interest in the dynamics of non-Abelian 
plasmas at very high temperatures or high densities. Due to the asymptotic freedom of 
quantum chromodynamics, one expects that quarks and gluons are no longer confined under 
such extreme conditions, but rather behave as free entities forming the so-called quark-gluon 
plasma. Within the next few years, a lot of efforts will be given to detect experimentally this 
new state of matter using heavy-ion colliders. Another domain of application concerns the 
physics of the early universe. If baryogenesis can finally be understood within an electroweak 
scenario, an understanding of the physics of the electroweak model in the high temperature 
regime where the spontaneous broken symmetry is restored, is essential for a computation 
of the rate of baryon number violation. 

It is therefore mandatory to devise reliable theoretical tools for a quantitative description 
of non-Abelian plasmas both in or out of equilibrium. While some progress has been achieved 
in the recent years [0,0, we are still far away from having a satisfactory understanding of 
the relevant relaxation and transport processes in non-Abelian plasmas, in particular when 
it comes to out-of-equilibrium situations. 

There are different approaches in the literature to study non-Abelian plasmas, ranging 
from thermal field theory to quantum transport equations or lattice studies. Even in the 
close-to-equilibrium plasma, and for small gauge coupling, the situation is complicated due 
to the non-perturbative character of long-wavelength excitations in the plasma. Most at- 
tempts to tackle this problem are based on a quantum field theoretical description of the 
non-Abelian interactions [Q. It has been conjectured that the plasma close to equilibrium 
allows for a description in terms of soft classical fields, as the occupation number for the 
soft excitations are large. While a classical transport theory approach for the non-Abelian 
case [^ is known to reproduce the one-loop thermal effective action [Q, it has never been 
exploited in full detail. The opposite holds true for Coulomb plasmas, where all the es- 
sential transport phenomena have been studied longly using techniques developed within 
(semi-) classical kinetic theory 0, while a quantum field theoretical approach has been un- 
dertaken only recently. 

Our approach aims at filling this gap in the literature of classical non-Abelian plas- 
mas. Here, we follow the philosophy of Klimontovitch [0], and our equations can be seen 
as the generalization of classical kinetic theory for Abelian plasmas to non-Abelian ones. 
Our essential contribution is considering the non-Abelian colour charges as dynamical vari- 
ables and introducing the concept of ensemble average to the non-Abelian kinetic equations. 
Equally important is the consistent treatment of the intrinsic non-linearities of non-Abelian 
gauge interactions. The observation that Klimontovitch's procedure leads to the Balescu- 
Lenard collision integral for Coulomb plasmas has motivated earlier derivations of similar 
(semi-) classical kinetic equations for non-Abelian plasmas [p|-p!0[]. However, these imple- 



mentations are not fully consistent, and have never been worked out in all generality. 

The starting point for a classical transport theory of non-Abelian plasmas is consider- 
ing an ensemble of classical point particles carrying a non-Abelian charge. They interact 



through self-consistent fields, that is, the fields generated by the particles themselves. The 
microscopic dynamics is governed by the classical equations of motion given by the Wong 
equations |[TT|]. When the number of particles is large, one has to abandon a microscopic 
description of the system in favour of a macroscopical one based on an ensemble average 
of all the microscopic quantities. This leads naturally to a description in terms of aver- 
aged quantities, and their statistical fiuctuations. By averaging the microscopic dynamical 
equations, we obtain effective transport equations for mean quantities. These contain the 
collision integrals of the macroscopic Boltzmann equation, which appear in this formula- 
tion as statistical correlators of fluctuating quantities. By subtracting the exact microscopic 
equations from the mean ones, we obtain the dynamical equations for the fluctuations them- 
selves. In principle, these two set of equations should be enough to consider all the transport 
phenomena in the plasma. 

This method is then applied in full detail to a thermal non-Abelian plasma close to equi- 
librium, which allows to employ some approximations. For a small plasma parameter, the 
two-particle correlators are small and can be neglected. In the case of small fluctuations, 
the dynamical equations simplify considerably. These conditions are always met for a small 
gauge coupling parameter, which, for simplicity, will be assumed throughout. Our approxi- 
mate equations are the leading order ones in a consistent expansion in the gauge coupling. 
However, we shall also see that the condition for a kinetic description to be valid could also 
be met for large gauge couplings. After taking statistical averages, we are able to explicitly 
integrate-out the fluctuations with momenta about the Debye mass. This gives the collision 
integrals which appear in the transport equations for the mean flelds. In addition to the 
dissipative processes in the plasma described by the collision term, we are able to deduce 
the stochastic source which prevents the system from abandoning equilibrium. This is an 
important result, because it allows to prove explicitly that the fluctuation-dissipation theo- 
rem holds, when switching from a microscopic to a macroscopic description of the system. 
These flndings for classical plasmas can be generalized to the case of quantum ones. The 
resulting dynamical equation match perfectly the effective theory for the ultra-soft modes 
as found by other approaches [T^^H • 



The lesson to be learned is thus two-fold: There exists a fully self-contained formalism to 
study classical non-Abelian plasmas in the flrst place, which opens in particular a door for 
applications to out-of-equilibrium situations. Second, this approach is -technically speaking- 
much easier than approaches based on the full quantum fleld theory. Some of the intrinsic 
complications of a quantum fleld theoretical description (like gauge-flxing, ghost degrees 
of freedom) can be avoided, and in the close to equilibrium plasma, the same effective 
dynamical equations are obtained. 

The paper is organized as follows. We begin with a review of the microscopic picture, 
based on the classical equations of motions for coloured point particles (sect. |I|). Chang- 
ing to a macroscopic description needs the introduction of a statistical average, which also 
allows the computation of correlators of fluctuations (sect. |T|). This procedure is then ap- 
plied to the flelds and the distribution function as to obtain dynamical equations for their 
mean values and their fluctuations. The dynamical equations are given in its most general 
form. Possible approximation schemes are detailed, and the interpretation of statistical cor- 
relators in terms of collision integrals is given (sect. fV\ ). The consistency of the procedure 



with the requirements of gauge invariance is shown for the complete set of equations, and 
some approximations to them (sect. ^. In order to apply the formalism to a plasma close 
to equilibrium, we discuss first the relevant physical scales for both classical and quantum 
plasmas (sect. |V| ). This is followed by a fully detailed derivation of the mean field dynam- 
ical equation for classical plasmas, which includes the integrating-out of the fluctuations 
with momenta about the Debye mass, and the computation of the collision integral and the 
related noise variable in a leading logarithmic approximation. A brief discussion of (itera- 
tive) solutions to the Boltzmann-Langevin equation for different momentum regimes is also 
presented (sect. [VI I[ ). We argue that these results can be translated to the case of quantum 



plasmas and detail the necessary changes. Some comments on related work are added as 
well (sect. |VIII|) . Finally, we present our conclusions (sect. |^, deferring to the appendices 
some technical details regarding the Darboux variables for SU(A^) colour charges (appendix 
^, and the derivation of a useful algebraic identity (appendix 0). 

II. MICROSCOPIC EQUATIONS FOR NON-ABELIAN CHARGED PARTICLES 

Let us consider a system of particles carrying a colour charge Q"^, where the colour index 
runs from a = 1 to A^^ — 1 for a SU(A^) gauge group. Within a microscopic description, 
the trajectories in phase space are known exactly. The trajectories x(r),p(r) and Q{t) for 
every particle are solutions of their classical equations of motions, the Wong equations |[TT[| 



m^=p^', (2.1a) 

dr 

m^ = gQ^'FrPu , (2.1b) 

(IT 

m^ = -gr"^p'^AlQ'^. (2.1c) 

Here, A^ denotes the microscopic gauge field. The corresponding microscopic field strength 
F" and the energy momentum tensor of the gauge fields 6^^^ are given by 

F;,[A] = d.A';, - d^Al + gr''^A\Al , (2.2) 

e^'^[A] = \g^''F;;^F^'' + f^t;; " (2.3) 

and f"''^'^ are the structure constants of SU(A^). We set c = fc^ = ^ = 1 and work in natural 
units, unless otherwise indicated. Note that the non-Abelian charges are also subject to 
dynamical evolution. Equation ( |2.1c| ) can be rewritten as D^-Q = 0, where Dj- = ^^D^ is 
the covariant derivative along the world line, and -D"^[A] = 8^5'^^ + gf°-'^'^A^ the covariant 
derivative in the adjoint representation. With Qa and F" transforming in the adjoint rep- 
resentation, the Wong equations can be shown to be invariant under gauge transformations. 
The equation ( p.lc|) ensures the conservation under dynamical evolution of the set of A^ — 1 
Casimir of the SU{N) group.f] 



•fFor SU{2), it is easy to verify explicitly the conservation of the quadratic Casimir QaQa- For 
SU{3), both the quadratic and cubic Casimir dabcQaQbQc, where dabc are the symmetric structure 



The colour current associated to each particle can be constructed once the solutions of 
the Wong equations are known. For every single particle it reads 

dr— Q,{r) 6'~'\x - x{t)) . (2.4) 



Employing the Wong equations ( |2.1|) we find that j'^ is covariantly conserved, Dfj_j^ = 
|TT| . Similarly, the energy momentum tensor associated to a single particle is given by [11 



dT^p''{T)6^^\x-x{r)) , (2.5) 

dr 

It is convenient to describe the ensemble of particles introducing a phase space density 
which depends on the whole set of coordinates x^,p^ and Qa- We define the function 

n{x,p,Q) = Y.jdT5^'\x-Ur)) S^'\p - Mr)) 5^'''-^\Q - Q,{r)) , (2.6) 

i 

where the index i labels the particles. This distribution function is constructed in such a 
way that the colour current 

J,^(x) =gjd^p rf(^'-i)Q ^Q, n(x, p, Q) (2.7) 

coincides with the sum over all currents associated to the individual particles, J^ = J2iJa, 
and is covariantly conserved, D^J'^ = 0. It is convenient to make the following changes in 
the choice of the distribution function. We will define a new function f{x,p,Q) such that 
the physical constraints like the on-mass shell condition, positive energy and conservation 
of the group Casimirs are factored out into the phase space measure. We introduce the 
momentum measure^ 

dP = d^p 2e{po)6{p^ - m^) , (2.8) 

and the measure for the colour charges 

dQ = d'Qcn5iQaQa-q2) , (2.9) 

in the case of SU(2). For SU{3) the measure is 

dQ = d^'Q CR6{QaQa - q2)S{datcQV'Q' " Qs) ■ (2.10) 



constants of the group, are conserved under the dynamical evolution. The last conservation can be 
checked using ( p.lc| ) and a Jacobi-like identity which involves the symmetric date c^^d antisymmetric 
fabc constants. 

^Note that in [|l|, we used a slightly different normalization of the measure. Here, the measure 
has an additional factor of (27r)^. 



For SU{N) , N — 1 (5-functions ensuring the conservation of the set of A^ — 1 Casimirs have to 
be introduced into the measure dQ. We have also introduced the represent at ion- dependent 
normahzation constant cr into the measure dQ, which is fixed by requiring J dQ = 1 (see 
appendix ^. The constant C2 is defined as 

dQQaQb = C26ab , (2.11) 

and depends on the group representation of the particles. 
With these conventions the colour current (12.71) reads now 



J!:{x)=gJdPdQp>^Q,f{x,p,Q) , (2.12) 

while the energy momentum tensor associated to the particles obtains as 

t'^^(x) = f dPdQ p^'p" f{x,p,Q) . (2.13) 

We will now come to the dynamical equation of the microscopic distribution functions 
n{x,p, Q) and f{x,p, Q), which will serve as the starting point for the subsequent formalism. 
The dynamical equation for n{x,p,Q) is the same as for f{x,p,Q). This is so because the 
physical constraints which we have factored out as to obtain f{x,p,Q) are not affected by 
the Wong equations. Employing ( ^.1[ ), we find 



which can be checked explicitly by direct inspection of (|2.6|) into (|2.14a|) (see appendix A of 



). In a self-consistent picture this equation is completed with the Yang-Mills equation, 

(D^F^aix) = J:{x) , (2.14b) 



and the current being given by ( p.l2| ). It is worth noticing that ( p.l4| ) is exact in the 



sense that no approximations have been made so far. The effects of collisions are included 
inasmuch as the Wong equations do account for them, although ( ^.14aD looks formally like 
a collisionless Boltzmann equation. 

For the energy momentum tensor of the gauge fields we find 

^^Q^"' = -F^/r^ . (2.15) 

On the other hand, using (|2.14a ) and the definition ( p.l3| ) we find that 

d,t^^ = F7J; . (2.16) 

which establishes that the total energy momentum tensor is conserved, d^{Q^'^ + 1^'^) = 0. 
To finish the review of the microscopic description of the system, let us recall the gauge 
symmetry properties of the set of equations ( |2.14|) (see for more details). From the 



definition of f{x,p, Q) it follows that it transforms as a scalar under a (finite) gauge trans- 
formation, f'{x,p,Q') = f{x,p,Q). This implies the gauge covariance of ( p.l4b| ) because 
the current ( |2.7D transforms like the vector Qa in the adjoint. The non-trivial dependence 
of / on the non-Abelian colour charges implies that the partial derivative 9^/ does not yet 
transform as a scalar, but rather its covariant derivative -D^/, which is given by 

D,[A]fix,p,Q) ^ [d,-gr'^QA,,,d^]fix,p,Q) . (2.17) 



Notice that ( |2 . 1 7| ) combines the first two terms of (|2.14a| ). (Here, and in the sequel we use 



the shorthand notation (9^ = d/dx'^, d^ = d/dp^ and 5^ = d/dQ"".) The invariance of 
the third term in (|2.14a| ) follows from the trivial observation that QaF'^^ is invariant under 



gauge transformations. This establishes the gauge invariance of (|2.14a|) . 

III. STATISTICAL AVERAGES 

If the system under study contains a large number of particles it is impossible to follow 
their individual trajectories. One has then to switch to a statistical description of the 
system. 

In this section, we describe in detail the statistical average to be used in the sequel. 
As we are studying classical point particles in phase space, the appropriate statistical av- 
erage corresponds to the Gibbs ensemble average for classical systems. We will review the 
main features of this procedure, defined in phase space. Let us remark that this deriva- 
tion is completely general, valid for any classical system, and does not require equilibrium 
situations. 

We will introduce two basic functions. The first one is the phase space density function 
M which gives, after integration over a phase space volume element, the number of particles 
contained in that volume. In a microscopic description it is just a deterministic quantity, 
and a function of the time t, the vectors x and p, and the set of canonical (Darboux) 
variables cf) and tt associated to the colour charges Q (see appendix ^). For SU{N), there 
are A^(A^ — l)/2 pairs of canonical variables, which we denote as = (0i, . . . ,(Pn(n-i)/2) 
and TT = (vTi, . . . , vrjv(Ar_i)/2)- The microscopic phase space density is given by 

Ar(x, p, 0, tt) = 5: 5(3) (x - Xi(t)) 6^'\p - Ut)) m - kit)) 5{tz - 7r,(t)) , (3.1) 

i 

where the sum runs over all particles of the system, and (xj, pj, 0^, ttj) refers to the trajectory 
of the i-th particle in phase space. Then ATrfx rfp dcpdn gives the number of particles at time 
t in an infinitesimal volume element of phase space around the point z = (x, p, <p, tt) 

Let us now define the distribution function JF of the microstates of a system of L identical 
classical particles. Due to Liouville's theorem, dj^/dt = 0. Thus, it can be normalized as 

/ dzidz2 . . . dzL T{zx, Z2, . . . , ZL,t) = 1 . (3.2) 

Here Zi denotes all the phase space variables associated to the particle i. For simplicity 
we have considered that there is only one species of particles in the system, although the 
generalization to several species of particles is rather straightforward. 

6 



The statistical average of any function Q defined in phase space is given by 

{Q) = / dzidz2 ...dzL G{zi, Z2,..., zl) T{zi, ^2, • • • , 2;^, t) ■ (3.3) 

The one-particle distribution function is obtained from JF as 

fi{zi,t) = V dZ2...dzLJ^{zi,Z2,...,ZL,t) , (3.4) 

where V denotes the phase space volume. Correspondingly, the two-particle distribution 
function is 

f2{zi, Z2, t) = V'^ dZ3... dZL J^{zi, Z2,..., Zl, t) , (3.5) 



and similarly for the /c-particle distribution functions. A complete knowledge of JF would 
allow us to obtain all the set of (/i, f2, ■ ■ ■ , fi) functions; this is, however, not necessary for 
our present purposes. 

Notice that we have allowed for an explicit dependence on the time t of the function 
JF, as this would be typically the case in out of equilibrium situations. We will drop this t 
dependence from now on to simplify the formulas. 

Using the above definition one can obtain the mean value of the microscopic phase space 
density. Microscopically N'{z) = J2i=i S{z — Zi), where Zi describes the trajectory in phase 
space of the particle i. The statistical average of this function is 

r ^ L 

{Af{z)) = / dzidz2 ■ . . dzLJ^{zi, Z2,..., Zl) Y1 S{z - Zi) = — fi{z) . (3.6) 

The second moment {Af{z)Af{z')) can equally be computed, and it is not difficult to see 
that gives 

W^m^')) = ^K^ - ^')/i(^) + ^^^^^h{z,z') . (3.7) 

Let us now define a deviation of the phase space density from its mean value 

5N{z) = M{z) - {M{z)) . (3.8) 

By definition {6Af{z)) = 0, although the second moment of this statistical fiuctuation does 
not vanish in general, since 

{6^{z)6^{z')) = {M{z)M{z')) - {M{z)){N{/)) . (3.9) 

If the number of particles is large, L ^ 1, we have 

{6U{z)6U{z')) = (^) 6{z - z')h{z) + {^)\2{z,z') , (3.10) 

where the function g2{z,z') = f2{z,z') — fi{z)fi{z') measures the two-particle correlations 
in the system. For a completely uncorrelated system g2 = 0. 



Notice that the above statistical averages are well defined in the thermodynamic limit, 
L,V —>■ oo but L/V remaining constant. Higher order moments and higher order correlators 
can as well be defined. 

We finally point out that the the function A/'(x, p, 0, tt) agrees with the microscopic 
function f{x,p,Q) introduced earlier, except for a representation-dependent normalization 
constant cr as introduced in sect. |I|. We will also swallow the density factors L/V into the 
mean functions /. Those small changes in the normalizations allow to simplify slightly the 
notations of the equations. 

IV. AVERAGING THE MICROSCOPIC EQUATIONS 

A. The mean fields and the fiuctuations 

In this section we perform the step from a microscopic to a macroscopic formulation of 
the problem. Using the prescription explained in sect. |T|, we take statistical averages of the 
microscopic equations ( p.l4|) . This implies that the distribution function f{x,p,Q), which 



in the microscopic picture is a deterministic quantity, has now a probabilistic nature and 
can be considered as a random function, given by its mean value and statistical (random) 
fiuctuation about it. Let us define the quantities 

f{x, p, Q) = fix, p, Q) + 6 fix, p, Q) , (4.1a) 

J,^(x) = J,^(x) + <5J„^(x), (4.1b) 

Alix) = Alix) + alix), (4.1c) 

where the quantities carrying a bar denote the mean values, e.g. / = (/), J = (J) and 
A = (A), while the mean value of the statistical fiuctuations vanish by definition, {6f) = 0, 
(6 J) = and (a) = 0. This separation into mean fields and statistical, random fluctuations 
corresponds effectively to a split into low frequency (long wavelength) modes associated to 
the mean fields, and high frequency (short wavelength) modes associated to the fiuctua- 
tions.Q We also split the field strength tensor as 

^;. = i^;. + /;. , (4.2a) 

^;. = F^ulA] , (4.2b) 

/;, = iD,a, - D,a,r + ^/""^« , (4.2c) 

using also D^ = Dn[A\. The term /" contains terms linear and quadratic in the fiuctuations. 
Note that the statistical average of the field strength (F^j^) is given by {F^^^) = F^^ + 
gf°'^'^{o!'^al), due to quadratic terms contained in f/^^. 

In the same light, we split the energy momentum tensor of the gauge fields into the part 
from the mean fields and the fiuctuations, according to 



**In the close-to-equilibrium plasma (sect. [VII ), we identify the relevant momentum scales 
explicitly. 
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Q,u ^ yf^-F^^Fr + K'F; ' , (4.3b) 

9'^^ = y^-'F^jr + Frc + kt,, + Wfp.^ar^" + r/i; " , (4.3c) 

The term 6^" contains the fluctuations up to quartic order. Due to the non-hnear character 
of the theory, we find that the ensemble average of the energy momentum tensor reads 

B. Dynamical equations for the mean fields and the fluctuations 

We perform now the step from the microscopic to the macroscopic Boltzmann equation, 
taking the statistical average of ( p. 141) . This yields the dynamical equation for the mean 



values, 

P" [d, - gQaP-^^d;) f = (r/) + (0 . (4.4a) 



We have made use of the covariant derivative of / as introduced in ( p.l7[ ). The macroscopic 
Yang-Mills equations are obtained as 

D.P"" + ( JL) = r . (4.4b) 

In ( [1.4|) , we collected all terms quadratic or cubic in the fluctuations into the functions 
ri{x,p,Q),^{x,p,Q) and Jfmcix). They read explicitly 

vix,p,Q)^gQaP^d;r^Jfix,p,Q) , (4.5a) 

ax,p,Q)^9p''r'V'{d^al 6fix,p,Q) +ga;ald;fix,p,Q)) , (4.5b) 

JZix) = 9 {f""'Dii,a,,X + r'^a,^, [{D^a'' - D^a^c + gr'^a^X^ } (x) . (4.5c) 



The corresponding equations for the fluctuations are obtained by subtracting ( [1.4D from 
( pA^ ). The result is 

p'^ [d^ - gQaF;;,d;) 6f = gQa{D,a, - D^a^Yp^^dlf 

+9P'a,,,r''Q,d^f + V + ^-{v + (4.6a) 

[D'a^ - D>^{D,a')y + 2g r'^^"" a,,, + J^ - (J^) = ^^"'^ • (4.6b) 

The above set of dynamical equations is enough to describe all transport phenomena in the 
plasma. 

While the dynamics of the mean fields ( |4.4| ) depends on correlators quadratic and cubic 
in the fluctuations, the dynamical equations for the fluctuations ([4.6|) do in addition depend 
on higher order terms (up to cubic order) in the fluctuations themselves. The dynamical 
equations for the higher order correlation functions are contained in (|4.6| ). To see this, 
consider for example the dynamical equation for the correlators {6f 6f). After multiplying 
( [4.6a| ) with 6f and taking the statistical average, we obtain 



P'^ {d, - gQaF^^d;) {6f6f) = gQalfdlf {{D^a, - D^a^fdf) 

+ gp'-r'^Qcd^f {a,,, 6f) + iiv + Sf) ■ (4.7) 

In the same way, we find for {6f 6f 6f) tlie dynamical equation 



P' [D, - gQaF;,d;) {6f 6f 6f) = gQaP'dlf {{D.a, - D^a.f 6f 5f) - {v + i^f Sf) 

+gp'r'Vcd'^f {ai>,, Sf 6f) + {{v + OSf Sf) , (4.8) 

and similarly for higher order correlators. Typically, the dynamical equations for correlators 
of n fluctuations will couple to correlators ranging from the order (n — 1) up to order (n + 2) 
in the fluctuations. From cubic order onwards, the back-coupling contains terms non-linear 
in the correlation functions.[^ 

The dynamical equation for the energy momentum tensor of the gauge fields obtains 
from the average of ( p.l5[ ). The corresponding one for the particles is found after integrating 



( p. 14a ) over dPdQp'^. They read 



a.e'^^ + d, {en = -Fr-ia - ifrsjua) - (/d Ju , (4.9) 

da''^ = F^Jua + ifrSJua) + {fn Jt , (4.10) 

such that the total energy momentum tensor is conserved. 

Finally, the condition for the microscopic current conservation translates, after averaging, 
into two equations, one for the mean fields, and another one for the fluctuation fields. From 
(D^J^) = we obtain 

iD,jna + gfabc{al6r^'') = 0. (4.11) 

For the fluctuation current, we learn from D^J^ — {D^j_J^) = that 

{D,6Jna + gfabc ialJ^ + al5Ji! - {a^dJ^)) = . (4.12) 



In sect. |V IJ| , it is shown that ( ^.llj ) and ( |4.12| ) are consistent with the corresponding equa- 



tions as obtained from the Yang-Mills equations. 

C. Second moment approximation and small coupling expansion 

The dynamical equations, as derived and presented here, are exact. No approximations 
have been performed. In order to solve them, it will be essential to apply some approx- 
imations, or to find a reasonable truncation for the hierarchy of dynamical equations for 
correlator functions. Here, we will indicate two approximation schemes, the second moment 



T^TThe resulting hierarchy of dynamical equations for the correlators is very similar to the BBGKY 
hierarchy. While the BBGKY hierarchy links the dynamical equations for n-particle distribution 
functions with each other, here, we rather have a hierarchy for the correlator functions. 
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approximation and the small coupling expansion. Although they have a distinct origin 
in the first place, we will see below (sect. ^ that they are intimately linked due to the 
requirements of gauge invariance. 

The so-called second moment approximation (sometimes also referred to as the polariza- 
tion approximation |^) is employed in order to simplify the dynamics of the fluctuations. It 
consists in equating 

V={V) , e = (0 , Jnu. = (Jfl.c) . (4.13) 

The essence of this approximation is that the dynamical equations for the correlators be- 
comes homogeneous. It is easy to see that ( [4.7|) or ( [4.8|) depend only on quadratic or cubic 



correlators, respectively, once ( [4.13 ) is imposed. This approximation allows to cut the infi- 



nite hierarchy of equations down to a closed system of differential equations for both mean 
quantities and statistical fluctuations. The mean fields couple to quadratic correlators, and 
all correlators of degree n couple amongst each others. This turns the dynamical equation 
for the fluctuations ([4.6|) into a differential equation linear in the fluctuations. This ap- 
proximation is viable if the fluctuations and the two-particle correlators are small (see also 
sect. ^. 

A seemingly different approximations concerns the non-Abelian sector of the theory. 
We shall perform a systematic perturbative expansion in powers of the gauge coupling g, 
keeping only the leading order terms. This can be done, because the differential operator 
appearing in the Boltzmann equation ( p. 14a ) or ([4.4a ) admits such an expansion. In a small 



coupling expansion, the term gQaFl^^d^ is suppressed by a power in g as compared to the 
leading order term p'^D^. Notice, that expanding the g appearing within p'^D^ of ( p.l7| ) is 
not allowed as it will break gauge invariance. In this spirit, we expand as well 

f = f^'^+9F'^+9'f^'^ + ... (4.14) 

and similarly for Sf. This is at the basis for a systematic organization of the dynamical 
equations in powers of g. To leading order, this concerns in particular the cubic correlators 
in (?7) and (Jhuc), which are suppressed by a power of g as compared to the quadratic ones. 
In principle, after these approximations are done, it should be possible to express the 
correlators of fluctuations appearing in ( [4.4| ) through known functions. This requires finding 
a solution of the fluctuation dynamics first. 



D. Collision integrals 

Let us comment on the interpretation of (rj) and (^) as collision integrals of the macro- 
scopic Boltzmann equation. The functions (77) and (^) appear only after the splitting ( |4.1| ) 
has been performed. This introduces new terms in the corresponding Boltzmann equa- 
tion ( |4.4| ) for the mean fields, which are interpreted as effective collision integrals for the 
macroscopic transport equation. In this formalism, the collision integrals appear naturally 
as correlators of statistical fluctuations. The fluctuations in the gauge fields cause random 
changes in the motion of the particles, while random changes in the distribution function of 
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the particles induce changes in the gauge fields. This is having the same effect as collisions, 
and yields a precise recipe as to obtain collision integrals from the microscopic theory. 
In this light, the second moment approximation ( [4.13[ ) can be interpreted as neglecting 



the back-coupling of collisions to the dynamics of the fluctuations. Also, neglecting cubic 
correlators appearing in {rj) or (Jane) in favour of quadratic ones, to leading order in an 
expansion in the gauge coupling and for small fluctuations, is interpreted as neglecting the 
three-particle collisions in favour of two-particle collisions. 

In order to flnd explicitly the corresponding collision integrals for the non-Abelian 
plasma, one has to solve flrst the dynamical equations for the fluctuations in the back- 
ground of the mean flelds. This step is interpreted as integrating-out the fluctuations. In 
general, this is a difficult task, in particular due to the non-linear terms present in ( [4.6|) . As 
argued above, this will only be possible when some approximations have been performed. 

In the Abelian limit, ( [4.4| ) and ( |4.6| ) reduce to the known set of kinetic equations for 
Abelian plasmas |^. In this limit, only the collision integral (77) survives. Here it is known 
that {f]) can be expressed explicitly as the Balescu-Lenard collision integral, after solving 
the dynamical equations for the fiuctuations and computing the correlators involved [0]. 
This proofs in a rigorous way the correspondence between fiuctuations and collisions in an 
Abelian plasma. 

V. CONSISTENCY WITH GAUGE SYMMETRY 

In this section we shall discuss the consistency of the present approach with the require- 
ments of the non-Abelian gauge symmetry. This discussion will concern the consistency 
of the general set of equations. The question of consistent approximations will be raised 
as well. In this section, we shall for convenience switch to a matrix notation, using the 
conventions A = AHa, Q = Q"'ta etc., as well as [ta,tft] = fahct^ and Tr tah = —\5ab- 

A. Gauge covariance of the macroscopic equations 

As a consequence of the Wong equations being gauge invariant, we already established 
in sect. ^ that the microscopic set of equations ( 2.14| ) transforms covariantly under (finite) 
gauge transformations 

gA^ = U{x){d, + gA,)U-\x) , U{x) = e^^-gt\x)ta , (5.1) 

with parameter e"(x). We also have f'{x,p,Q') = f{x,p,Q), and 

Q' = U{x)QU-\x), d'Q = U-\x)dQUix) , F'^^ = U{x)F,,U-\x) . (5.2) 

The question raises as to which extend this symmetry is conserved under the statistical 
average, performed when switching to a macroscopic description. This concerns in particular 
the subsequent split of the gauge field into a mean (or background) field, and a fiuctuation 
field 

A^ = A^ + a^. (5.3) 
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This separation is very similar to what is done in the background field method (BFM) [p!5| 
Two symmetries are left after the splitting is performed, the background gauge symmetry, 



gA'^ = U{x){d, + gA,)U-\x) , a'^ = U{x) a^U'^x) 
and the fluctuation gauge symmetry, 

gA'^ = 0, ga'^ = U{x) (d, + g{A^ + a^)) U'^x) . 



(5.4) 



(5.5) 



Under the background gauge symmetry, the fiuctuation field transforms covariantly (as a 
vector in the adjoint). In the first step, we will split ( |2.14a| ) according to ( |5.3| ). It follows 
trivially, that the resulting equation is invariant under both ( |5.4| ) and ( p.5|) , if both / and 
5f transform as /, that is, as a scalar. 

The next step involves the statistical average. We require that the statistical average of 
the fiuctuation vanishes, (a) = 0. This constraint is fully compatible with the background 



gauge symmetry, as (a) = is invariant under ( |5.4| ). Any inhomogeneous transformation 
law for a, and in particular ( ^.5|) , can no longer be a symmetry of the macroscopic equations 
as the constraint (a) = is not invariant.]^ 

We rewrite now the transport equations in matrix convention. We have 



p^ [d, + 2g Tr {QF,,) d;) f = (v) + {0 ■ 



(5.6a) 
(5.6b) 



with 



v{x,p,Q) = -2g TTiQf,,)p^d;6f{x,p,Q) , (5.7a) 

e(x,p,Q) = -2gp^^ TT{[Q,d'^]a^)6f{x,p,Q) -2/ Tr {[a^,a,]Q) p^^d;f{x,p,Q), (5.7b) 

JL(x) = g [5^ K, a"]] + [a^, [D>^, a^] - [D^ a'^]] + g'[a^, K, a'']]. (5.7c) 

and 



p^ [D, -gTi (QF,,) d;) 6f = -2gTT iQ[D„ a,] - Q[D,, a^^p^dlf 

-2(?p^Tr([g,9«]a^)/ + r/ + e-(r/ + , (5.8a) 

D,, [D^, a^]] - \d^, [D,, all + MF'^, «.] + JL " (JL) = ^J' ■ (5.8b) 



It is now straightforward to realize that (|5.6| ) - ( ^.8|) transform covariantly under (|5.4| ). It 
suffices to employ the cyclicity of the trace, and to note that a^ and background covariant 
derivatives of it transform covariantly. 



WThis is similar to what happens in the BFM, where the fluctuation gauge symmetry can no 
longer be seen once the expectation value of the fluctuation field is set to zero. However, the 
symmetry ( ^.5| ) will be observed in both ( [4.4a ) and ( [4.6a| ), as long as the terms linear in (a) are 
retained. 
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B. Consistent current conservation 



In ( 4.11|) and ( 4.12| ), we have given the equations which imply the covariant current 
conservation of the mean and the fluctuation current. However, this information is contained 
both in the transport and in the Yang-Mills equation. It remains to be shown that these 
equations are self-consistent. 

We start with the mean current J . Performing g J dPdQ Q of the transport equation 
( p.6a| ), we find 



= [D^,J^]+g{[a„6J^]) . 
This is ([4.11|) . Here, we used that 

ldPr]{x,p,Q) = 0, 
JdPp'^F,,d;f{x,p,Q) = 0, 

j dPdQQ i{x,p,Q) = -\a^,5J^] 
g fdPdQQp^^DJ{x,p,Q) = [D^,J^] . 



Taking the background covariant derivative of ( |5.6bD , we find 



(5.9) 

(5.10a) 
(5.10b) 
(5.10c) 
(5.10d) 

(5.11) 



which has to be consistent with ( ^.9|) . Thus, combining these two equations we end up with 
the consistency condition 



0=[D„{J[:J]+g{[a„SJ^]) . 
In the appendix P, we established the identity 

= [D,, JL] + 9K, 6J^] + g[a,, (J^)] 



(5.12) 



(5.13) 



which follows using the explicit expressions for 6J'^ from ( ^.8b|) and for J^^^ from ( ^.7c|) . Tak- 
ing the average of ( |5.13D reduces it to ( [5. 12] ), and establishes the self-consistent conservation 
of the mean current. 

The analogous consistency equation for the fluctuation current obtains from ( |5.8a| ) after 
performing g J dPdQ Q, and reads 



= [D^, 5J^] + g[a^, 6J^^] + g[a,, J^] - g{[a^, 6J^]) 
This is ( ^.121 ). Here, in addition to ( |5.10| ), we made use of 

2gj dPdQQ Tr ([g,9«] a^) f{x,p,Q) = g[a^, J''] . 

The background covariant derivative of (5.8b) is given as 



(5.14) 



(5.15) 
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= [D„ 5J^] + g [a., [D„ Fn\ - [D„ J^] + [5„ ( JL)] • 
Subtracting these equations yields the consistency condition 

= [D,, JL] + qK, 5 J'] - [D,, (JL)] - 9{K, ^r]) + qK, J'] - g fa., [D,, Fn 



(5.16) 



(5.17) 



Using ( ^.6b|) , (|5.12| ) and (|5.13| ) we confirm ( |5.17| ) exphcitly. This estabhshes the self- 
consistent conservation of the fiuctuation current. 



C. Gauge-consistent approximations 

We close this section with a comment on the gauge consistency of approximate solutions. 
The consistent current conservation can no longer be taken for granted when it comes to 
finding approximate solutions of the equations. On the other hand, finding an explicit 
solution will require some type of approximations to be performed. The relevant question 
in this context is to know which approximations will be consistent with gauge invariance. 

Consistency with gauge invariance requires that approximations have to be consistent 
with the background gauge symmetry. From the general discussion above we can already 
conclude that dropping any of the explicitly written terms in ( [4.4|) - ([4.6| ) is consistent 
with the background gauge symmetry ( [5.41 ). This holds in particular for the polarization or 
second moment approximation ( [4.13| ). 

Consistency of the second moment approximation ( [4.13 ) with covariant current conser- 
vation turns out to be more restrictive. Employing 7^^,^ = {J^^^) implies that ( [5.12[ ) is only 
satisfied, if in addition 







Du, (K 



/l; [^ ,0- 



(5.1^ 



holds true. This is in accordance with neglecting cubic correlators for the collision integrals. 
Similarly, the consistent conservation of the fiuctuation current implies the consistency 
condition ([5.17[) , and holds if 



= [a„ (JL)] 



(5.19) 



vanishes. It is interesting to note that the consistent current conservation relates the second 
moment approximation with the neglection of correlators of gauge field fiuctuations. We 
conclude, that ( [4.13[ ) with ( [5.18[ ) and ( p.l9[ ) form a gauge-consistent set of approximations. 



VI. PHYSICAL SCALES: CLASSICAL VERSUS QUANTUM PLASMAS 

In the remaining part of the paper, we study the classical and quantum non-Abelian 
plasma close to equilibrium. Prior to this, we shall present a discussion of the relevant 
physical scales of both relativistic classical and quantum plasmas, close to equilibrium. We 
will restore here the fundamental constants h, c and fc^ in the formulas. 
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To discuss the relevant physical scales in the classical non-Abelian plasma, it is conve- 
nient to discuss first the simpler Abelian case, which has been considered in detail in the 
literature ^. At equilibrium the classical distribution function is given by the relativistic 
Maxwell distribution, 

ripo) = Ae-P"/'^^^ , (6.1) 

where A is a dimensionful constant which is fixed once the mean density of particles in the 
system is known. For massless particles po = pc, and if no further internal degrees of freedom 
are present, the mean density is given by iV = 87TA{kBT/cY. The interparticle distance is 
then f ~ iV~^/^. As we are considering a classical plasma, we are assuming f ^ AdB, where 
AdB is the de Broglie wave length, AdB ~ ^/p, with p some typical momenta associated to 
the particles, thus p ~ ksT /c. The previous inequality implies therefore A <^ 1/h , which 
is the condition under which quantum statistical effects can be neglected. 

There is another typical scale in a plasma close to equilibrium, which is the Debye length 
r^j. The Debye length is the distance over which the screening effects of the electric fields 
in the plasma are felt. For an electromagnetic plasma, the Debye length squared is given 

rl = ksT/A'KNe'^ . (6.2) 

Notice that the electric charge contained in the above formula is a dimensionful parameter: 
it is just the electric charge of the point particles of the system. 

In the classical case, and in the absence of the fundamental constant h, the only dimen- 
sionless quantity that can be constructed from the basic scales of the problem is the plasma 
parameter e. The plasma parameter is defined as the ratio 

e = fVr|) . (6.3) 

The quantity 1/e gives the number of particles contained in a sphere of radius r^. If e <^ 1 
this implies that a large number of particles are in that sphere, and thus a large number 
of particles are interacting in this volume, and the collective character of their interactions 
in the plasma can not be neglected. For the kinetic description to make sense, e has to 
be small [0. This does not require, in general, that the interactions have to be weak and 
treated perturbatively. 

Let us now consider the non-Abelian plasma. The interparticle distance is defined as 
in the previous case. The main differences with respect to the Abelian case concerns the 
Debye length, defined as the distance over which the screening effects of the non-Abelian 
electric fields in the plasma are noticed. It reads 



D 



kBT/A7rNg^C2 , (6.4) 



where C2, defined in ( [^.11| ), is a dimensionful quantity, carrying the same dimensions as the 
electric charge squared in ( |6.2|) . The coupling constant g is a dimensionless parameter. In 
the non-Abelian plasma one can also construct the plasma parameter, defined as in 
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It is interesting to note that there are two natural dimensionless parameters in the non- 
Abehan plasma: e and g. The condition for the plasma parameter being small translates 
into 

(!f)"^"v«i, (M 

which is certainly satisfied for small gauge coupling constant g <^ 1. But it can also be 
fulfilled for a rarefied plasma. Thus, one may have a small plasma parameter without having 
a small gauge coupling constant. This is an interesting observation, since the inequalities 
e <^ 1 and (? -C 1 have different physical meanings. A small gauge coupling constant allows 
for treating the non-Abelian interaction perturbatively, while e -C 1 just means having a 
collective field description of the physics occurring in the plasma. In principle, these two 
situations are different. If we knew how to treat the non-Abelian interactions exactly, we 
could also have a kinetic description of the classical non-Abelian plasmas without requiring 

Now we consider the quantum non-Abelian plasma, and consider the quantum coun- 
terparts of all the above quantities, as derived from quantum field theory. It is com- 
mon to change the phase space measure d^xd^p to the standard quantum normalization 
d^xd^p/{27rh)^, which makes the quantum distribution function dimensionless.^^ For a quan- 
tum plasma at equilibrium the one particle distribution function is 

^b'(Po) = gp„/,,T _ 1 ' /f"(Po) = ^,,/,,r ^ ^ , (6.6) 

where the subscript B/F refers to the bosonic/fermionic statistics. For a plasma of massless 
particles the mean density is iV ~ {ksT /TicY . The interparticle distance f ~ N~^/^ becomes 
of the same order as the de Broglie wavelength, which is why quantum statistical effects can 
not be neglected in this case. 

The value of the Debye mass is obtained from quantum field theory. It depends on the 
specific quantum statistics of the particles and their representation of SU{N). From the 
quantum Debye mass one can deduce the value of the Debye length, which is of order 

(6.7) 

It is not difficult to check that the plasma parameter, defined as in ( |6.3| ), becomes propor- 
tional to g^. Thus e is small if and only if g ^ 1. This is so, because in a quantum field 
theoretical formulation one does not have the freedom to fix the mean density N in an ar- 
bitrary way, as in the classical case. This explains why the kinetic description of a quantum 
non-Abelian plasma is deeply linked to the small gauge coupling regime of the theory. 




^=As done in the standard textbooks the factor (27r?i)^ is introduced into the measure, although 
some authors introduce it into the distribution function. We follow the first option. 
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VII. THE CLASSICAL PLASMA CLOSE TO EQUILIBRIUM 

In this section we put the method to work for a hot non-Abehan plasma close to equi- 
librium. A prerequisite for a kinetic description to be viable is a small plasma parameter 
e <^ 1. We shall ensure this by imposing a small gauge coupling constant f? -C 1. Then, all 
further approximations as detailed in the sequel can be seen as a systematic expansion in 
powers of the gauge coupling. 

A. Non-Abelian Vlasov equations 

We begin with the set of mean field equations ( [4. 41 ) and neglect the effect of statistical 



fluctuations entirely, 6f = 0. In that case, ( [4.4|) becomes the non-Abelian Vlasov equations 



p^[D^-gQaF^,d';)f = 0, (7.1a) 

D^F^"" = r , (7.1b) 

where the colour current is given by 

J'a{^)=9 E [dPdQQaP''f{x,p,Q), (7.1c) 



hclicitics 
species 



We will omit the species and helicity indices on the distribution functions, and in the sequel, 
we will also omit the above sum, in order to keep the notation as simple as possible. Equation 
( [7.1a| ) is then solved perturbatively, as it admits a consistent expansion in powers oi g. Close 
to equilibrium, we expand the distribution function as in ( |4.14| ) up to leading order in the 
coupling constant 

f{x,p,Q) = r{po)+9pHx,p,Q) . (7.2) 

In the strictly classical approach, the relativistic Maxwell distribution ( p.l|) at equilibrium 
is used. Here, we consider only massless particles, or massive particles with ?7i -C T, such 
that the masses can be neglected in a first approximation. We will consider also internal 
degrees of freedom, two helicities associated to every particle. 

It is convenient to re-write the equations in terms of current densities. Consider the 
current densities 

J'a,-aA^,P) = 9 P' I dQQa,---QaJ{x,p,Q), (7.3a) 

Ja',-a„ix,v) = I dP JI...,SX,P) . (7.3b) 



Here, v^ = (l,v) with v^ = 1. The measure dP integrates over the radial components. It 
is related to ( p.8| ) by dP = dPdQ/An, and reads 

dP = Att dpo d\p\ IpP 2e(po) 5(/) (7.4) 



for massless particles. The colour current is obtained performing the remaining angle inte- 
gration J{x) = J ^J{x^ v). From now on we will omit the arguments of the current density 
JT", unless necessary to avoid confusion. 

We now insert (|7.2|) into (|7. 1|) and expand in powers of g. The leading order term 
p ■ Df^^'{po) vanishes. After multiplying (4.4a) by gQaP''/po, summing over two helicities, 



and integrating over dPdQ, we obtain for the mean current density at order g 

v^D^J" + mlv^v^F^o = , (7.5a) 

D^F^'' = r , (7.5b) 



with the Debye mass 



m 



,dn{p) 
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STrg'C, r dpp'^LM , (7.6) 

JO dp 



The solution to ( [7.5a|) is now constructed with the knowledge of the retarded Green's func- 
tion 

tv^'D^G,,t{x,y;v)=6^'^\x-y) . (7.7) 

It reads 

Grctix, y; v)ab = -i9{xo - yo)5^^'^ (x - y - v(xo - yo)) Uabix, y) , (7.8) 

where IJab is the parallel transporter obeying v^D^^ Uab{x, y)\y=x-vt = 0, and Uab{x, x) = 6ab- 
One finds 

-^ / dr Uabix, X - VT)v^v''F,o,bix - vt) . (7.9) 

An JO 

The above colour current agrees with the hard thermal loop (HTL) colour current @,§], 
except for the value of the Debye mass. 

From ( [7.5a ) it is easy to estimate the typical momentum scale of the mean fields. If 



the effects of statistical fluctuations are neglected (and as we will see, this is equivalent to 
neglecting collisions), the typical momentum scales associated to the mean current and the 
mean field strength are of the order of the Debye mass m£,. We will refer to those scales 
as soft scales. The momentum scales with momenta ^ itld will be referred to as ultra-soft 
from now on. 



B. Leading order dynamics for mean fields and fluctuations 

We now allow for small statistical fluctuations 6f{x,p,Q) around ( |7.2| ), writing 

f{x,p,Q) = ripo)+gP\x,p,Q)+6f{x,p,Q) (7.10) 



and re-write the approximations to ( [4.4|) and ( |4.6|) in terms of current densities and their 
fluctuations. Note that the fluctuations 6f{x,p,Q) in the close to equilibrium case are 
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already of the order of g. This observation is important for the consistent approximation in 
powers of the gauge couphng. As a consequence, the term gf^^^ in (|7.1CI| ) will now account 



for the ultra-soft modes for momenta <^ mr). Integrating-out the fluctuations results in an 
effective theory for the latter. 

As before, we obtain the dynamical equation for the mean current density at leading 
order in g, after multiplying ( [4.4a| ) by gQaP'^/po, summing over two helicities, and integrating 
over dPdQ. The result is 

v^-D^JP + mlv'^v^F^o = iv') + ie) , (7.11a) 

D.F'^'^ + iJL) =r . (7.11b) 

In a systematic expansion in gf, we have to neglect cubic correlator terms as compared to 
quadratic ones, as they are suppressed explicitly by an additional power in g. Therefore, we 
find to leading order 



r/,^ = g fdP^iD^a, - D,a,f d; SJ^,{x,p) , (7.12a) 
•^ Po 

ea = -9fabcv^a'^6j'''', (7.12b) 

JPf^ = gf^^(^Dfa^aP + 5^'^al(D^'aP - DPa^'Y) . (7.12c) 



The same philosophy is applied to the dynamical equations for the fluctuations. To leading 
order in g, the result reads 

[v^D^ 6JP)^ = -mlv^v^ [D^ao - D^a^) ^ - gfabcV^alJ^^'^ , (7.13a) 

V'' {d^dac^M + gA^ {fa^, 6m + fbmdSac)) ^J^d = 9^^0.7 Umac ^M + fmM^ac) J^d , (^-ISb) 

DV-D^pa)) +2gUcFrac,u = SJ^ . (7.13c) 



The typical momentum scale associated to the fluctuations can be estimated from (|7.13| ). 



We find that it is of the order of the Debye mass ~ m^), that is, of the same order as 
the mean fields in (|7.5| ). This confirms explicitly the discussion made above. The typical 
momentum scales associated to the mean fields in ( 7.11 ) is therefore <^ mi). 



C. Integrating-out the fluctuations 



We solve the equations for the fluctuations (|7.13|) with an initial boundary condition for 



6f, and a^(t = 0) = 0. Exact solutions to ( [7.13aD and ( |7.13bD can be obtained. 
Let us start by solving the homogeneous differential equation 

v^D^6JP = 0, (7.14) 

with the initial condition 6J'/^{t = 0,x, f). It is not difficult to check, by direct inspection, 
that the solution to the homogeneous problem is 

SJP{x, v) = Uab{x, X - vt) 5Ji;{t = 0, X - vt, v) . (7.15) 

20 



The solution of ( [7.13a] ) is now constructed using the retarded Green's function (|778|). For 
Xo = t > the complete solution can be expressed as 



5JI^{x, v) = - j dr Uab{x, Xr) i mlv'^v^ (^D^ao - Doa^) (xr) + gfMcV^a'l{xr)Jl^{xr, v) 

+Uai,{x,xt)5j(;{xuv). (7.16) 

We have introduced Xr = x — vt, and thus Xt = (0,x — vt). Since a^(t = 0) = 0, one can 
check that the above current obeys the correct initial condition. 

The equation ( [7.13b| ) can be solved in a similar way. The solution is 



SJ^f,{x, V) = Uam{x, Xt)Ubn{x, Xt) 6J^^{Xt, v) 

dr Uam{x, Xr)Ubn{x, Xr) {fmpc^nd + fnpd^mc) V^aF{Xr)Jcd{Xr, v) . (7.17) 



Now we seek for solutions to the equation ( [7.13c| ) with the colour current of the fluctuation 



as found above. However, notice that this equation is non-local in a^, which makes it 
difficult to find exact solutions. Nevertheless, one can solve the equation in an iterative 
way, by making a double expansion in both gA and gj- This is possible since the parallel 
transporter U admits an expansion in gA, so that the current 5 J'' can be expressed as a 
power series in gA 

5JP = SJP^""^ + 5JP^^'^ + 5JP^^^ + • • • , (7.18) 



and thus ( |7.13c| ) can be solved for every order in gA. To lowest order in gA, using Uab 



^ab + 0{gA), equation (|7.13c| ) becomes 



d' [d.a^l - d.a^Sa) = ^4'! ■ (7-19) 

Using the one-sided Fourier transform 0, and ( [7.16| ), we find 



6J!ii'\k) = Kbik)a^>) - gfabc [^-^ I -^yaf\q) ^-(k - g, .) 



47r -i k-v J (27r)4 p 

dQ, 6JP{t = 0,k,v) ^ .720) 

Att —ik-v 

where 11^^ (/c) is the polarization tensor in the plasma, which reads 

Kbik) = S^bml l^-gP^ + ^0/^^) , (7.21) 

and agrees with the HTL polarization tensor of QCD 0,^, except in the value of the Debye 

mass. Retarded boundary conditions are assumed above, with the prescription ko — >• kQ+iO~^. 

We solve ( [7.19| ) iteratively in momentum space for a^ as an infinite power series in gS, 

aT = «r°^ + «r'^ + <°"^ + • • • (7-22) 
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where the second index counts the powers of the background current gj. Notice that in 
this type of Abehanized approximation, the equation ( [7.19|) has a (perturbative) Abehan 



gauge symmetry associated to the fluctuation a^. This symmetry is only broken by the 
term proportional to J in the current. It is an exact symmetry for the term ajj^'^-* in the 
above expansion. We will use this perturbative gauge symmetry in order to simplify the 
computations, and finally check that the results of the approximate collision integrals do 
not depend on the choice of the fluctuation gauge. 

Using the one-sided Fourier transform, we find the following results for the longitudinal 
fields, in the gauge k ■ a^°'°'' = 0, 

while we find 

—k^ + 11t J 4:71 —ik-v 



for the transverse fields.Q The functions lii/Tik) are the longitudinal/transverse polariza- 



tion tensor of the plasma, -PjT(k) = 5ij — kikj/'k'^ the transverse projector, and af = P^aj. 
In the approximation g <^ 1, it will be enough to consider the solution of leading (zeroth) 
order in gA, and the zeroth and first order in gJ'. The remaining terms are subleading in 
the leading logarithmic approximation. However, notice that we have all the tools necessary 
to compute the complete (perturbative) series. If we could solve equation ( |7.13c|) exactly, 
it would not be necessary to use this perturbative expansion. 



D. The statistical correlator of fluctuations 



With the explicit expressions obtained in (|7.20|) , (|7.23|) and (|7.24|) , we can express all 



fluctuations in terms of initial conditions 6J'l^{t = 0,x, f) and the mean fields. 

In order to compute the correlator of initial conditions, we will make use of the result 
obtained in sect. ^. For each species of particles or internal degree of freedom, the statistical 
average over initial conditions can be expressed as 

{6 fit = 0, X, p, Q) 6 fit = 0, x', p\ Q')) = 6^'^ (x - x')^^^) (p - p')5(Q - Q')fix, p, Q) 

+^2(x,p,Q;x',p',Q'), (7.25) 



***In [ffl], we used a more condensed notation. There, the functions aj ^ and a\ ^ have been 
denoted a\ ^_ and a- :J_. 
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where the function g2 obtains from the two-particle correlator, and 

6{Q - Q') = —6{(j) - ci>') 6{tt - 7t) , (7.26) 

Cr 

and 0, TT are the Darboux variables associated to the colour charges Qa- The appearance 
of the factor 1/cr in the above expression is due to the change of normalization factors 
associated to the functions M and /, as we mentioned at the end of sect. |ITl| . The above 
statistical average is all we need to evaluate the collision integrals below. 



From ( [7.25| ) one deduces the statistical average over colour current densities 6J. We 



expand the momentum (5-function in polar coordinates 

5(3)(p - pO = 4 Kp - P') 5^'Hfiv - fiv) , (7.27) 

where Q^ represents the angular variables associated to the vector v = p/|p|. After simple 
integrations we arrive at 



+ ~g^l,{^,v;^',v'), (7.2^ 



Be = IQn' / dpp' r^^ip) . (7.29) 

Jo 

The function ^2^ is obtained from the two-particle correlation function ^2- Notice that we 
have neglected the piece gf^^^ above, as this is subleading in an expansion in g. 

Since we know the dynamical evolution of all fluctuations we can also deduce the dy- 
namical evolution of the correlators of fluctuations, with the initial condition ( |7.25| ). This 



[6j;{t = 0, X, v) 5Jt{t = 0, x', v')) = 2g'Bc C^ 5"^ v^ 5(3)(x - x') 6^'\n^ - l^v 
where v'^ = (l,v), and 



corresponds to solving (4/7) in the present approximation. It is convenient to proceed as 



follows 0. We separate the colour current ( |7.16| ) into a source part and an induced part 



5J^ = (5:^1..,, + (5 j;::,^, . (7.30) 

The induced piece SJ'induccd is the part of the current which contains the dependence on a^, 
and thus takes the polarization effects of the plasma into account. The source piece SJ''^ 

A -C^ IT *~^ source 

is the part of the current which depends only on the initial condition, given by the solution 
of the homogeneous equation ( [7.14| ). This splitting will be useful since ultimately all the 
relevant correlators can be expressed in terms of correlators of SJ'J^^^^^. 

From the explicit solution ( |7.1(j| ) and the average ( [7. 28] ) we then find, at leading order 
in g and neglecting the non-local term in ( |7.28| ) 



SJ:SJ^, v) 5Jt:.Jx', v')) = 2g'Bc C2V^v"'S^'\^ - x' - v(t - t'))5^'\n^ - fi^O 

xU'^'ix, X - vt)U^\x\ x' - v't') . (7.31) 

Here, and from now on, we neglect the non-local piece gf'^jj- It can be shown []^ that they 
give contributions which decrease rapidly with time, so that for asymptotic large times, they 
vanish. 
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Expanding the parallel transporter U, and switching to momentum space we find the 
spectral density to zeroth order in gA 

(SJ^ SJX^''^]^'^ = '^a'Bc C25''\vl5^^\n^ - n^,)i2n)Sik ■ v) . (7.32) 

As an illustrative example, let us compute the correlator of two transverse fields a. Using 
( [7.24a| ) and (|7.32|) one arrives at 

aT\k) air\q)) = g'Bc C,5.,{2.Y5^^\k + g)^^^ J^v.v, 5{k ■ v) . (7.33) 

Since the imaginary part of the polarization tensor, which describes Landau damping, can 
be expressed as 

ImliZi^) = -Satmlrrko J^v'^v-' 6{k ■ v) , (7.34) 

the statistical correlator can finally be written as 

(C'°'(*) -■r(l)) = '-^^^^f^i2''f^''''ik + ,) . (7.35) 

Here, we have used the relation 

2g^C2Bc = ATrTml . (7.36) 



Equation ( |7.35| ) is a form of the fiuctuation dissipation theorem (FDT), which links the 



dissipative processes occurring in the plasma with statistical fiuctuations. 



E. The collision integral 

We are now ready to compute at leading order in g the collision integrals appearing on 
the r.h.s. of (|7.11a|) . We shall combine the expansions introduced earlier to expand the 
collision integrals in powers of jf (while retaining only the zeroth order in gA), 

(^) = (^(o)) + (^(i)) + (^(2)) + ... , (7.37) 

and similarly for {rj) and {Jruc)- We find that the induced current {JLI) vanishes, as do 
the fiuctuation integrals (r/^^^) and {^^^^). The vanishing of {JaJc) is deduced trivially from 
the fact that (a^°'°^a|, ) ^ 5ab, while this correlator always appears contracted with the 
antisymmetric constants fabc in <^uc- To check that {r]^'^'>) = 0, one needs the statistical 
correlator {SJ^^SJ^^), which is proportional to Y^a^aab = for SU{N). The vanishing of 
(^(0)) jg consistent with the fact that in the Abelian limit the counterpart of {r]) vanishes at 
equilibrium [^. Finally, {^^^^) = due to a contraction of fabc with a correlator symmetric 
in the colour indices. 
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In the same spirit we evaluate the terms in the colhsion integrals containing one J field 
and no background gauge A fields. Consider 



+ 9fccieV''l°^dTj^4xr,v) (a^ (x) a'^^f (Xr)) }• (7.38) 

Using the values for a^ and 5J'^^' as found earlier, we obtain in momentum space 

(ej:i(^,^)) ^ -g'C^NBcvJ^Civy) {J^{k,v)-J^{ky)) , (7.39) 



with 

C(v,vO 






i2ny 



ViP[jiq)v'j 



-q" + Ht 



{2-K)5{q ■ v){2Tx)5{q ■ v') . (7.40) 



Here, the symbol ~ means that only the leading terms have been retained. To arrive at the 
above expression we have used the SU{N) relation fabcfabd = Ndd- Within the momentum 
integral, we have neglected in the momenta of the mean fields, k, in front of the momenta of 
the fiuctuations, q. As we discussed above, the momenta associated to the background fields 
are much smaller than those associated to the fiuctuations which justifies this approximation 
to leading order. This is precisely what makes the collision integral, which in principle 
contains a convolution over momenta, local in A;-space (resp. x-space). The only remaining 
non- locality stems from the angle convolution of ( [7.4U| ). Notice that we have only written 
the part arising from the transverse fields a, as the one associated to the longitudinal modes 
is subleading. This is easy to see once one realizes that the above integral is logarithmic 
divergent in the infrared (IR) region, while the longitudinal contribution is finite. At this 
point, we can also note that the collision integral computed this way is independent on the 
(perturbative) Abelian gauge used to solve equation ( |7.19| ). This is so because the collision 
integral computed this way can always be expressed in terms of the imaginary parts of the 
polarization tensors ( |7.34| ) in the plasma, which are known to be gauge- independent. 

The integral (|7.40|) has also been obtained in |jl3|, on the basis of a phenomenological 



derivation of the Boltzmann collision integral for a quantum plasma. The only difference 
consists in the value of the Debye mass appearing in the polarization tensor. 

In any case, the transverse polarization tensor 11^ vanishes at go = 0, and the dynamical 
screening is not enough to make (|7.40|) finite. An IR cutoff must be introduced by hand 
in order to evaluate the integral. With a cutoff of order gmo we thus find at logarithmic 
accuracy 

<^(v.v')«^ln(l/.)-^il^ (T.41) 



Using also the relation ( [7.36| ) we finally arrive at the collision integral to leading logarithmic 
accuracy, 
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e«(x,.;)) = -^NT\n{l/g)v,J^I{v,v')X%x,v'), (7.42) 

J(v, v') = 5(2) (v - v') - /C(v, v') , /C(v, v') = 1 (l^lL , (7.43) 



where we have introduced 6^^\v - V) = A7i6^'^\Q^ -K), I l^^^^^v - V) = 1. 

We can verify exphcitly that the colhsion integral to leading logarithmic accuracy is con- 
sistent with gauge invariance. This should be so, as the approximations employed have been 
shown in sect. [VC| on general grounds to be consistent with gauge invariance. Evaluating 



the correlator in ( [4. Ill ) in the leading logarithmic approximation yields 

gUc{al{x)6Ji:{x)) = -^NTln{l/g) J ^^I{^.y)J^{x,v') , (7.44) 



which vanishes, because 



/^J(v,v')=0. (7.45) 

We thus establish that D^J^^ = 0, in accordance with ( [7.11b|) in the present approximation. 



F. The source for stochastic noise 

The collision integral obtained above describes a dissipative process in the plasma, so 
in principle, it could trigger the system to abandon equilibrium. Whenever dissipative 
processes are encountered, it is important to identify as well the stochastic source related to 
it. This is the essence of the fluctuation-dissipation theorem (FDT). Phenomenologically, 
this is well known, and sometimes used the other way around: imposing the FDT allows to 
add by hand a source for stochastic noise with the strength of its self-correlator fixed by the 
dissipative processes. 

In the present formalism, we are able to identify directly the source for stochastic noise 
which prevents the system from abandoning equilibrium. This proves, that the FDT does 
hold (analogous considerations have been presented in ||12[)- The relevant noise term is given 
by the contributions from the transversal gauge fields in i^^^^^ While its average vanishes, 
(^^(0)^ = 0, its correlator 

^a'-'Kx, v) C,^'\y, v')) = g'fa^jMeV^v"^ (<(a;) SJ:j:Jx, v) at{y) ^J-^d/, v')f^ (7.46) 

does not. In order to evaluate this correlator we switch to Fourier space. Within the second 
moment approximation we expand the correlator {Sf6f6f6f) into products of second order 
correlators {Sf6f){6f6f) and find 



Ca^'\k, V) e,^'\p, V')) = g'fa,jMeV^v'^ 



d'^q f d^r 



{2'kYJ (27r)4 
+ {a^^q) Sj}lTip - r, ^0) {SJ.^l'lfik - g, v) a^ir))} (7.47) 
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In the leading logarithmic approximation we retain only the contributions from the trans- 
verse modes. Evaluating the correlators leads to 

{^^^i^,v)e^^{y,v')) = ^ y ,^ ln{l/g)v^v"'I{^^y)6'^'6'^'\x-y) . (7.48) 

After averaging over the angles of v and v', and using the relation ( [7. 361 ), the correlator 
becomes 

^^j{x) e^^M = 2T "^ ^ NT \n{l/ 9)6'^' 6^^ 6^'Hx-y). (7.49) 



In particular, all correlators {^?o)i^)^{o)iy)) vanish. Equation ( [r.49| ) identifies Co){x) as a 
source of white noise. 

The presence of this noise term does not interfere with the covariant current conservation 
confirmed at the end of the previous section. This can be seen as follows. The noise 
term enters ( [7.44[ ) as the angle average over the 0-component of C,a)\{x,v). As we have 



established above, the logarithmically enhanced contribution from the noise source stems 
from its correlator (|7.48|) . Averaging the temporal component of (|7.48|) over the angles of 



V, and using ( 7.45|) , it follows that 

(e?o?(^,-')/^e(>,-)) = 0. (7.50) 

We thus conclude, that the temporal component of the noise, ^?qJx,v), has no preferred 
v-direction, which entails that / ^^PqJx, w) = in the leading logarithmic approximation. 
Thus, the mean current conservation is not affected by the noise term. 



G. Mean Field Equations and Non-Abelian Ohm's law 

We have managed to obtain the following set of mean field equations, after integrating- 
out the statistical fluctuations (from now on, we drop the bar to denote the mean fields) 



v^'D^J'ix, v) = -mlvPy^F^oix) - ^ v" f —^ J(v, V)J%x, v') + C{x, v) , (7.51a 

J An 

Air 



D,F>^'= l-^J'{x,v). (7.51b) 



Here, we denote by C^p{x,v) the stochastic noise term identified in the preceding section, its 
correlator given by ( [7.481 ). We also introduced 

7 = £:iVT In (1/^) , (7.52) 

which will be identified as (twice) the damping rate for the ultra-soft currents. We shall 
refer to ( [7.51a| ) as a Boltzmann-Langevin equation, as it accounts for particle interactions 



via a collision integral as well as for the stochastic character of the underlying fluctuations. 
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The Boltzmann-Langevin equation ( |7.51a| ) has three distinct scale parameters, the tem- 



perature T, the Debye mass ttid, and the damping term 7. In the leading logarithmic 
approximation, these scales are well separated, 7 -C rra/j <^ T. This is why ( [7.51[ ) is 
dominated by different terms, depending on the momentum range considered. For hard 
momenta, ( 7.51a| ) is only dominated by the l.h.s., reducing it to the (trivial) current of hard 



particles moving on world lines. For momenta about the Debye mass, the term proportional 
to m\) becomes equally important, while the noise term and the collision integral remain 
suppressed by '')/mD- The resulting current is then given by ( |7.9| ), the HTL current. Mo- 
mentum modes below the Debye mass are affected by the damping term of the collision 
integral. Close to the scale of the Debye mass, the higher order corrections to ( [7.91 ) are 
obtained as an expansion in 7/f ■ D, 

00 
J^{x,v) = Y.JUx,v), (7.53) 



n=0 



where the current densities J(t\{x^v) obey the differential equations 

v^D^Ji;^^{x, v) = -mlv'^v'^F^oix) + C{x, v) (7.54a) 

v^D,Ji;,^ix,v) = -7^'^/^X(v, v')J-(Li)(a:,t;') (7.54b) 



Apart from the noise term, the leading order term in this expansion, J'!'qJx,v), coincides 
with the HTL current ( |7.9| ). All higher order terms J'f^Jx,v) are smaller by powers of 
~ (7/f ■ -D)", and recursively given by 



J'^Q-^{x,v) = / dTU{x,x — VT)l—Tn'j^v''v^Ffj_o{x — vT)+(^''{x — vT,v)\ , (7.55a) 

J^^){x,v) = -7 I cirf/(x,x -H^'^y ^ J(v, v')J(Li)(x,t;') (7.55b) 

This expansion is consistent with covariant current conservation. For every partial sum up 
to order n, we have 



D, (^S) + Ji:i) + ■■■ + J^n)) = (7.56) 

This expansion describes correctly how the presence of the collision integral modifies the 
ultra-soft current. However, it has a limited domain of validity. In particular, the over- 
damped regime where w • D <C 7 can not be reached, because the effective expansion pa- 
rameter diverges. 

Alternatively, one can separate the local from the non-local part of the collision integral 
to perform an expansion in the latter only. The effective expansion parameter is then 
7/(f ■ D + 7), which has a better IR behaviour. We expect this expansion to be much 
better for the spatial than for the temporal component of J''{x,v). This is so, because the 
term proportional to the non-local part /C(v,v') of the collision integral in ( [7.51a]) gives no 



contribution to the dynamical equations of the spatial component J*, after angle averaging 
( [r.51a[ ) over the directions of v. However, for the dynamical equation of the temporal 
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component, this term precisely cancels the local damping term, which is of course a direct 
consequence of current conservation. 

In this light, we decompose the current as in ( [7.53D to find the differential equations 

{v^'D^ + 7) J'(^o)(^' ^) = -mlv'^v^F^oix) + ('{x, v) , (7.57a) 

{v'^D, + ^)Jf^^{x, v)=iv^j^ /C(v, v'):r(^i)(x, v'). (7.57b) 

The retarded Green's function of the differential operator 

tiv^'D^ + ^f) G,,,{x,y;v) = 6^^\x-y) , (7.58) 

reads, for t = Xq — z/q, 

Gretix, y, v)ab = -i9{t)5^'^^ (x - y - vt) exp(-7t) Uab{x, v) . (7.59) 

and the iterative solution to the Boltzmann-Langevin equation obtains a; 

j7(Q-)(a;, f) = / dr e'x.Y>{—'^T)U{x,x — vt) \—ra\,v''v^ Fjq{x — vt) + C,''{x — VT,v)\ (7.60a) 

<J 

/•oo p (jCi f 

JCn)i^,v)=^J^ dTexp{-^T)U{x,x-VT)v'J -^ IC{v,V)JI^^_,^{x - VT, v') . (7.60b) 

This expansion in consistent with current conservation, if the angle average of J'P^^ vanishes 
for some n. This follows from taking the temporal component of ( |7.57|) and averaging the 
equation over v, to find 

DoJL + D,Jio) = -iJL , (T.61a) 



for the individual contributions, and 



DoJl) + D,Ji. = 7Jf„-i) - iJin) , (7.61b) 



D, (^ro) + ^fi) + • • • + •/("„)) + lJ(n) = , (7.61c) 



for their sum, which is consistent if jJ^n) is vanishing for some n 



in) 

To leading order, the ultra-soft colour current JT/qJx, v) in ( |7.60a|) has the same functional 



dependence on the field strength and on the parallel transporter as the soft colour current 
( [7.9|) . There is, however, an additional damping factor e~'^'^ in the integrand. 

We shall now consider the overdamped regime (or quasi-local limit) of the above equa- 
tions. Consider the mean field currents ( |7.60| ). The terms contributing to these currents 
are exponentially suppressed for times r much larger than the characteristic time scale I/7. 
On the other hand, the fields occurring in the integrand vary typically very slowly, that is 
on time scales <^ l/mD- Thus, the quasi-local limit consists in the approximation 



^I^^We thank D. Bodeker for pointing out a mistake in the reasoning presented in an earlier version. 
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Uab{x, X -vt) ^ Uab{x, x) = 6ab , Fjo{x - vt) ^ Fjo{x) . (7.62) 

In this case the remaining integration can be performed. The solution for the spatial current 
J*(x) stems entirely from the leading order term ( |7.60aD . All higher order corrections vanish, 
because they are proportional to 

/^vX:(v,v') = 0. (7.63) 

J 4:71 

Therefore, we obtain for the spatial current from ( [7.60| ) 

where a denotes the colour conductivity of the plasma. The noise term becomes 

K^) = -J^Cix,v) , (ulix) vi{y)) = 2TaS'^SabS^'\x-y) (7.64b) 

The noise term appearing in the Yang-Mills equation becomes white noise within this last 
approximation. The fluctuation-dissipation theorem is fulfilled because the strength of the 
noise-noise correlator ( [7.64b|) is precisely given by the dissipative term of (|7.64a|) . This is 



the simplest form of the FDT. The colour conductivity in the quasi-local limit has been 
discussed by several authors in the literature |[T3|JT( 



For the complete set of gauge field equations in the quasi-local hmit we still need to know 
J°(x). The iterative solution ( [7.60| ) gives J?n){^) = to any finite order. Therefore, we use 
instead the un-approximated dynamical equation for J'^{x,v), which yields, averaged over 
the directions of v, current conservation. In combination with the solution (|7.64a|) for the 
spatial current, the Boltzmann-Langevin equation ( |7.51| ) finally becomes 

D^F''' = aE' + v' , (7.65a) 

DiE' = ,f (7.65b) 

Dojo = -aJ° - Div' . (7.65c) 

It is worth pointing out that already in the leading logarithmic approximation the noise 
term appearing in the Yang Mills equation is not white, except in the local limit ( [7.64b| ). 
The noise in the Boltzmann-Langevin equation, on the other hand, is white (see ( [7.49|) ), 
when averaged over the directions of v. 

For numerical computations, which can in principle take into account the non-localities 
of the problem, it might be more convenient to work with the two set of equations (|7.51| ), 
rather than with a non-local stochastic gauge field equation. 

VIII. THE QUANTUM PLASMA CLOSE TO EQUILIBRIUM 
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A. The quantum plasma from transport theory 

Up to now we have made an entirely classical derivation of a Boltzmann equation with 
collision integrals and stochastic sources, and we have finally derived the mean gauge field 
equations. The basic ingredients for such a derivation were the classical equations of motion 
and the classical statistical averages introduced in sect. fT|. The following natural step is 
quantizing the whole procedure in order to obtain quantum Boltzmann equations and the 
corresponding mean gauge field equations. 

In order to quantize this formulation one has to abandon the concept of classical trajec- 
tories, and introduce commutators for all the canonical conjugate pairs of variables. One 
should also take into account that quantum particles are indistinguishable. The natural 
formulation of the quantum problem is then given in terms of Wigner functions and den- 
sity matrices. We will not present a rigorous discussion of the quantum counterpart of 
our formulation, but rather present a minimal set of changes in our equations which allow 
to consider also quantum plasmas close to equilibrium. We leave for a future project a 
much more rigorous discussion based on first principles of the quantum formulation of the 
problem. 

Our starting observation is that even in quantum plasmas the physics occurring at soft 
and ultra-soft scales can be encoded into classical or semiclassical equations. The reason 
for this is that the occupation number for soft modes close to equilibrium is very large, 
suggesting that a description in terms of classical equations might also be vahd to describe 
the physics of the soft and ultra-soft scales in the quantum plasmas. 

Therefore, in order to consider a quantum plasma, we shall need to make several changes. 
The first step consists in expanding the mean distribution function around the appropriate 
quantum statistical distribution function. For a plasma close to equilibrium, these are given 
by (|0|). We also change the normalisation of the distribution functions as indicated in 
sect. ^ 

This suffices to obtain the correct quantum value for the Debye mass, appearing in the 
Vlasov equation, and thus to reproduce fully the HTL effective theory in the leading order 
in^ 0. 

When fluctuations are considered as well, it is equally important to modify the classical 
correlator ( [7.25| ) to the corresponding quantum statistical one. For bosons, and for every 



internal degree of freedom, one has 

(5/x,p,Q 5/x'y,Q') = (27r)35(3)(x - ^')6^'\p - p')6{Q - g')/B(l + /b) 

+ ^2^(x,p,Q;x',p',g') , (8.1) 

while the corresponding correlator for fermions is 

(5/x,p,Q <5/x'y,Q') = (27r)35(3)(x - x')(5(3)(p - p')5iQ - Q')Ml - /p) 

+ ^2^(x,p,Q;x',p',Q') . (8.2) 

The functions ^2 ^ire related to the bosonic/fermionic two-particle correlation function. 
The above relations should be derived from first principles in a similar way as our equation 
(7.25). In the limit /b/f <^ 1 they reduce to the correct classical value. It has also to 
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be pointed out that the correlators ( |8.1| ) and ( |8.2| ) have been derived for the case of both 
an ideal gas of bosons and ideal gas of fermions close to equilibrium, matching the change 
described above. This can be taken as the correct answer in the case that the non-Abelian 
interactions can be treated perturbatively. 

With the above in mind, we can now describe the minimal set of changes to our com- 



putations of sect. [VI1| which allows to treat the quark-gluon plasma close to equilibrium. 
We will consider gluons in the adjoint representation with C2 = N, and Np quarks and 
Np antiquarks in the fundamental representation, with C2 = 1/2. All particles carry two 
helicities. We will neglect the masses of the quarks, m <^T. 
The value of the the quantum Debye mass squared becomes 

ml = -4 r dpp' (n^ + Np^] . (8.3) 



TT^ Jo \ dp dp 

Evaluating explicitly the integral, one finds ml = g'^T'^{2N + Np)/6. 

The correlator of colour currents densities are then modified, according to the changes 



mentioned above. Since we now consider different species of particles in (|7.28|) 6*2-80 is 
replaced by a sum over different species of particles. For the quark gluon plasma, all our 
equations of sect. |VI1| remain valid if we replace C2BC by 



2A^ f^ . o^„., ^„„. 2A^ 



C2Bc = — dpp'fl\l + Q) + '-^ dpp'fP{l-fP). (8.4) 

TT JO 71 Jo 



It is curious that the relation ( 7.36 ) for the quantum values of the above quantities remains 
unchanged, that is, 

2(^2 J2 C2BC = AnTml (8.5) 

species 

holds true also in the quark-gluon plasma. Since this combination appears in front of all 
our collision integrals, we find a universal value for the coefficient of ( |7.42| ), 7 for both 
the classical and quantum plasmas. This is always the case in the leading logarithmic 
approximation, if the IR cutoff used is of order gmp), where rrip) would correspond to the 
classical or quantum Debye mass, respectively. The value 7/2 can be identified with the 
damping rate of a hard transverse gluon Jl^ . 

With these observations, one does not need to repeat the computations that we per- 
formed in sect. |VII| . In particular, the final mean field equations of sect. |VII G| only change 



in the value of the Debye mass. 

B. Comparison to related work 

Let us briefiy comment on some related work. A similar philosophy to ours has already 
been followed by Selikhov |Q . He used the semiclassical limit of quantum transport equations 
for the Wigner functions associated to gluons and quarks, which reduce to our starting 
classical transport equations. He used a procedure of splitting both the Wigner functions 
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and vector gauge fields into mean values and statistical fluctuations. A key point is how 
the statistical correlator of fluctuations in a quantum framework can be derived. Selikhov 
relied on the same type of statistical correlator as derived in ( [7.281 ). However, it should be 
stressed that this statistical correlator is only correct in the pure classical framework, for 
classical statistics. This can not reproduce the correct prefactors of the quantum collision 
integral. Also, the FDT is not satisfied in this case. Instead, the correct correlators are given 
by (|8.1|) and ( ^.21) . Also, the colour current he found is not covariantly conserved. This is 
so because the non-local term in the collision integral ( [7.421 ), proportional to /C(v,v'), has 
been neglected. 

The first to derive the mean field equations (|7.51|) and the related noise correlator ( [7.48| ) 



for the quantum plasma was Bodeker [|12[. His approach uses the local version of the HTL 
effective action as the starting point to integrate-out the modes with momenta about the 
Debye mass. He profits from the observation that the soft field modes behave classically. 
This allows the definition of a classical thermal average with a weight given by the HTL 
Hamiltonian. Although the starting point and the techniques involved appear to be quite 
different from our approach, it remains intriguing that the integrating-out of momentum 
modes within a field theory results in the same effective theory as integrating-out statistical 
fluctuations within the present transport theory approach. 



A different line has been followed by Arnold, Son and Yaffe |13|, who realized that 
Bodeker's effective theory has a physical interpretation in terms of kinetic equations. They 
managed to derive the relevant collision term of the Boltzmann equation on phenomenolog- 
ical grounds. 

Very recently, the quantum collision integral has been obtained within a quantum fleld 



theoretical setting by Blaizot and lancu |14]. The derivation relies on a gauge covariant 



derivative expansion of the quantum fleld equations. While conceptually very different to 
our approach |^ , the approximations used in [|1^ are very similar to the ones we performed. 



This is maybe not too surprising after all, as the explicit computations in both approaches 
are based on a consistent expansion in powers of the gauge coupling constant. It seems only 
that the concept of statistical fluctuations has not been introduced in [Q , which may be a 
reason for why the source of stochastic noise, necessary for a correct macroscopic description 
of the plasma, has not yet been identifled. The stochastic noise can probably be derived by 
considering the effects of higher order correlation functions [TE\. 



Based on our approach |]1| , the quantum collision integral for the transport equation has 
been obtained as well by Valle |]TP|]. He started from the HTL effective theory and found the 



correct damping rate after imposing a fluctuation-dissipation relation. The noise term in 
his flnal equations is however missing, which again would entail that the system abandons 
equilibrium. 

Finally, it is interesting that the collision integral can be interpreted in terms of Feynman 
diagrams |12|JT^ . Bodeker also made a diagrammatic derivation of his effective theory [^ 



This is a much more lengthy and cumbersome task, and shows, on the other hand, the very 
efficiency of a kinetic approach, as it corresponds to a re-organisation of the perturbative 
series. 
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IX. DISCUSSION 

We have presented a self-consistent approach to study classical non-Abelian plasmas. 
Let us summarize here again our starting assumptions to derive the effective transport 
equations. A system of point particles carrying non-Abelian charges is considered. Their 
microscopic equations of motion are the Wong equations. In order to describe an ensemble 
of these particles, we introduced an ensemble average, which takes also the colour charges 
as dynamical variables into account. This yields finally a set of transport equations for 
both mean quantities and statistical fluctuations, and gives a recipe to obtain explicitly the 
collision integrals for macroscopic transport equations. This approach is consistent with 
gauge invariance, and admits systematic approximations. Most particularly, it does not rely 
on close-to-equilibrium situations. These techniques, applied since long to Abelian plasmas, 
have never been fully exploited for the non-Abelian case. Our approach is aimed at closing 
this gap in the literature of non-Abelian plasmas. 

We applied this method to non-Abelian plasmas close to thermal equilibrium. A suf- 
ficiently small gauge coupling parameter is at the basis for a systematic expansion of the 
dynamical equations. Neglecting fluctuations yields to leading order the known non-local 
expression for the soft current in terms of the soft gauge fields (HTL approximation). 
Integrating-out, in addition, the fluctuations to leading logarithmic order, that is with mo- 
menta about m£), results in a Boltzmann-Langevin equation for the ultra-soft modes, 

v^^D^J^ix^v) = -mlv^v'^F^oix) - 7 v^ J ^I{^r,V)j\x,v') + C{x,v) . 

It contains a collision term and a related noise term, with 7 = g'^NT In {l/g)/47r, while the 
stochastic source ( obeys 

{C^ix,v)Cay^v')) = ^lTmlv^v'^Iiv,v')6atS^'\x-y) . 

The above is completed with the Yang-Mills equations. Solving the Boltzmann-Langevin 
equation (see sect. [VII G|) , one obtains the Yang- Mills equation for the ultra-soft fields 



Surprisingly, the dynamical equations are the same for both classical and quantum plasmas, 
the only difference being the value for the Debye mass. This conclusion relies also on the use 
of an infrared cut-off of order grriD, where mo is the classical or quantum Debye mass. For 
the quantum plasmas, our result agrees with the quantum collision integrals found in the 
literature using different methods |T^|-Q . The main effect of the fluctuations with momenta 
about the Debye mass is the introduction of a damping term and a source of stochastic noise 
into the above expression. Note also that the damping coefficient 7 is the same for classical 
and quantum plasmas. 

Our work establishes a link even beyond the one-loop level between the classical transport 
theory approach as presented here, and a full quantum field theoretical treatment. It would 
be most desirable if this connection could be further substantiated. This should also yield 
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a quantitative criterion for the applicability of the -technically speaking- much simpler 
approach based on the classical point particle picture. 

Let us finally emphasize that the same IR problems, which are due to the unscreened 
magnetic modes, appear both for classical and quantum plasmas. This suggests that the 
solution for these IR divergences might also be the same in the two cases. Therefore, it 
seems profitable to seek for a solution to this problem in the much simpler framework of 
classical transport theory, rather than in a quantum field theoretical approach. 
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APPENDIX A: DARBOUX VARIABLES 



The statistical averages defined in sect. Ill have to be performed in phase space. The 



colour charges Qa are not real phase space variables 0. It is possible to define the set of 
Darboux variables associated to the Qa charges. 

For SU{2) we define the new set of variables 0, it, J by the transformation ^ 



Qi = cos (f)\\P^^ Q2 = smcpVJ"^ - TT^ Qa = TT (Al) 

where tt is bounded by —J < vr < J. The variables (p, vr form a canonically conjugate 
pair, while J is fixed by the value of the quadratic Casimir, which is constant under the 
dynamical evolution. One can define Poisson brackets with these canonical variables, under 
which the colour charges form a representation of SU{2), {Qa,Qb}pB = fabcQc- With the 
above change of variables, one can easily fix the value of the representation normalization 
constant cr introduced in (|2.9|) . From the condition ^ dQ = 1 one finds cji = 1/211^/0^. 
From the condition J dQQaQb = C2Sab one gets q2 = 36*2. This entirely fixes the value of 
cr as a function of C2. 

The Darboux variables associated to SU{3) were defined in 0, and will not be discussed 
explicitly here. 

We should also comment that in the pure classical framework, C2 carries dimensions of 
he. After quantization, the quadratic Casimirs should take quantized values proportional 
to h. The Poisson brackets then have to be replaced by commutators. 

APPENDIX B: CONSISTENT CURRENT CONSERVATION 

In this appendix we verify explicitly the identity 

= [D,, JU + 9K, Sr] + g[a,, ( J^)] , (Bl) 

which is at the basis for the proof of the consistent current conservation of both the mean 



field and the fluctuation current in sect. [VB| . The following check is algebraic, and it 
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will make use of symmetry arguments like the antisymmetry of the commutator and the 
tensors F^u,ffiu, andof the cychc identity [4, [tb,tc\] + [U, [tc,ta\] + [tc, [ta,tb]] = 0. The 
identity [D^, D^] = gF^^ is employed as well. To simplify the computation, we will seperate 
the fluctuation part of the field strength ( |4.2c| ) into the term linear and quadratic in a, 
according to 



J^J.v Jl,iJ.i' ~r J'. 



I2,^u ; 



fl,^^u = [Di,, ay\ - \D, 



Recall furthermore, using (|5.7cD and (|5.8b|) , that 



V') ^flj J 



V^l 



f2,^iu — g[0'^l■,a'u\ ■ 






(B2) 



(B3) 
(B4) 



are functions of the fluctuation field a. The first term of ( PTD reads, after inserting Jf^c from 



(B5) 



[D,, JL] = [Du, [D„ in] + 9[Du, K, in] + 9[D., K, fH] ■ 
Using SJ from (P^), it follows for the second term of ( PT| ) 

g[a„ SJ'^] = g'{a,, [a„ F'^l] + g[a,, [D„ /r]] + g K, JL] - 9 W. {JIJ] (B6) 

The last term of ( |B6D will be canceled by the last term in (|BI|). We show now that the first 
three terms of (p5|) and ( p6|) do cancel one by one. The first term in ( pS]) can be re-written 
as 



p/ifi 



P., [D,. /r]] = [[D., i^,], /n - [D., [D,, /ri] = w.,, /r 

Similarly, the first term of (p^ ) yields 
For the second term in (p5|) we have 



(B7) 



(B8) 



g\Dy, [a^, /(""]] = 5f[a^, {D^, /H] + gWv^ a^]' /i 



eii.v^ 



i/l^i-i 



-g[a^,[D^Jl 
which equals (minus) the second term of (p^). Finally, consider the third term of (|B6|) , 

9 [au, JL] = 9>u, [D„ K, a1]] + 9>u, K, r]] 



(B9) 



= -9[D^,,[a^J2'']] 
which equals (minus) the third term of (B5|). This establishes (|BT|). 
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